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Abstract 

In this paper, we study the dual representation for generalized multiple stopping problems, 
hence the pricing problem of general multiple exercise options. We derive a dual repre- 
sentation which allows for cashflows which are subject to volume constraints modeled by 
integer valued adapted processes and refraction perio ds modeled by st o pping times. A s 
such, this extends the works b y ISchoenmakersI (I20~10h . Bender (1201 la ). Bender ( 2011rJ ). 
Aleksandrov and Hamblv ( 201fj() . and Meinshausen and HamblvT l 2004 ) on multiple exercise 
options, which either take into consideration a refraction period or volume constraints, but 
not both simultaneously. We also allow more flexible cashflow structures than the additive 
structure in the above references. For example some exponential utility problems are covered 
by our setting. We supplement the theoretical results with an explicit Monte Carlo algorithm 
for constructing confidence intervals for the price of multiple exercise options and exemplify 
it by a numerical study on the pricing of a swing option in an electricity market. 



2010 AMS subject classifications: 60G40, 65C05, 91B25. 
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1 Introduction 

The last decades have seen ground breaking developments of Monte Carlo methods for American 
options based on multidimensiona l underlying price processes. In t he lat e nineties the regression 
based methods by Carriere ( 19961 ). Longstaff and Schwart2 ( 200ll ). and Tsitsiklis and Van Rov 
(|200ll ) may be considered as main breakthroughs. In general these methods provide lower bounds 
on the option price by constructing an approximation to the optimal exercise (stopping) time 
via regression on a set of basis functions . As such these appr oaches are termed " prima l" . At the 
beginning of this century Rogers! ( 2002 ) and independently Haugh and Koeanl ( 20041 ) provided 
the next breakthrough by presenting a "dual" representation for the optimal stopping problem 
corresponding to the American option pricing problem. In this representation the option price 
is expressed as infimum of an expecta tion over a set of ma r tinga les. (The key behind this dual 
representation can already be found in Davis and Karatzas ( 19941 ). in fact.) While in the primal 
methods the central problem is to find a "good" stopping time, in the dual problem one needs to 
find a "good" martingale, which leads to an upper bound for the price of an American option. As 
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one of the standard numeri cal approaches to compute du al upper bounds for American options 
by Monte Carlo we refer to lAndersen and Broadid (|2004l ). 

During the same time, in the emerging electricity markets products with a multiple of exercise 
opportunities, such as "swing options", became popular. Naturally, pricing of such a product 
leads to a multiple stopping problem, and so numerical methods for solving multiple stopping 
problems were called for. In this respect, generalization of the existin g primal regression methods 
for sta ndard optimal stopping was just a matter of routine. Further. iBender and Schoenmakers 
developed a kind of policy iteration for multiple stopping. However, regarding the dual 
approach the situation was not so clear. Meinshausen and Hamblv ( 20041 ) proposed a dual 
representation for the multiple stopping problem via expressing the excess value due to each 
additional exercise right by an infimum of an expectation over a set of martingales and a set 
of sto ppin g times. This l ine of research was carried out further by lAleksandrov and Hamblv 
(|2010h and lBenderl d2011bll in the context of dual pricing of multi-exercise options under volume 
constraints. Recently, Schoenmakers ( 2010l ) introduced a dual representation for the price of 
a multiple exercise option in contrast to the dual representation for the excess value of an 
additional right. This new dual representation involves an innmum over martingales only and 
can thus be considered as a more natural extension of the dual representation for single exercise 
options. This approach was generalized by Bender ( 2011al ) to a continuous time setting involving 
(constant) refract ion periods. 

In the meantime iKobvlanski et al.l (|201ll ) introduced and studied multiple stopping problems 
in the primal sense in a far more general context, where the payoff is considered to be some 
abstract functional of an (ordered) sequence of stopping times. The goal of the present paper is 
to find (pure) martingale dual representations for such generalized multiple stopping problems 
in a discrete time setting. As we will show, such representations can be constructed even in 
a most general setting. However, for practical implementation these general representations 
unfold their full strength only, if applied to some more specifically structured cashflows. In this 
respect we study a generic payoff structure with both multiplicative and additive structure that 
incorporates (integer valued) volume constraints and refraction periods given by stopping times. 
We furthermore provide an explicit Monte Carlo based algorithm and give a detailed numerical 
study exemplifying the pricing of swing options. Comparing to existing works, the numerical 
experiments reveal that, by and la rge, the dual algorithms due to ou r ne w representation s applied 
to the problem type considered in Aleksandrov and Hamblv ( 2O10l ) and Bender ( 2011b ) produce 
tighter upper bounds on the option price, in particular when the number of exercise rights is 
large. We moreover present a numerical example which involves swing options subject to both 
volume constraints and refraction periods and give tight confidence intervals for the respective 
option prices. We underline that the latter example cannot be treated by the dual methods 
presented in the literature so far. 

The structure of the paper is as follows: In Section 2 we derive a dual representation for general 
multiple stopping problems in terms of a family of martingales. As this family of martingales 
is typically too large for practical purposes in general, we specialize to a generic cashflow with 
additive and multiplicative structure which incorporates volume constraints and refraction pe- 
riods in Section 3. In this section we then prove two dual representations. One is in terms of 
the Doob decomposition of the Snell envelopes of an auxiliary family of stopping problems, the 
other one only requires approximations of these Snell envelopes. In Section 4 we explain how to 
build a Monte Carlo algorithm for computing confidence intervals on the value of the multiple 
stopping problems based on the results of Section 3 and perform some numerical experiments 
in the context of swing option pricing. 
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2 General multiple stopping problem 



In this section we consider a multiple stopping problem in discrete time i = 0, . . . , T, where T € N 
is a fixed and finite time horizon. We further introduce a "cemetery time" d := T + 1 where all 
rights will be exercised, which are not exercised up to time T. For a given filtration (J-i)o<i<d and 
a number L of exercise dates we next consider a cashflow X as a map 
which satisfies for all < i\ < • • • < il < d, 

Xi 1) ... ) i L is Ti L -measurable, 
E l^n,...,^! < oo. 

Now consider the stopping problem (sup:=ess.sup, Ej := Ejr.), 

Y* L = sup E 4 X rli ... iTi , 

where the supremum runs over a family of ordered stopping times T k ,1 < k < L. 
Let us define for k = 2, . . . , L and < ji . . . < jk-i < v < d, 

y> : ' H ' - -P VrX h _ jk _ 1>rK ^, (1) 

r<T k <---<T L 

with the convention that for k = 1, we put l^*^' := and for = L + l, we put l^.* ' 51 '"' = 

■*-3\,...,j L - 



Proposition 1. We have the following reduction principle 



Y, 



*L-k+l,j 1 ,...,j k _ 1 



sup E r y T 



*L-k,ji,...,j k ^ 1 ,T 



r > jk-1- 



(2) 



Proof. This principle can be straightforwardly proved in an inductive manner, but it can also 
be considered as a disc rete time version of a related result in a continuous time setting from 
Kobvlanski et al.l (|201lh . □ 



In what follows the following remark turns out to be useful. 

Remark 2. We say that a martingale (M r ) r > p is a Doob martingale of (Y r ) r > p , whenever there 
exists a predictable process (A r ) r > p , such that Y r — M r + A r is T v -measurable for any r > p. In 
particular, for any two Doob martingales (M r ) r > p and (M r ) r > p of (Y r ) r > p , it holds 

r-l 

M r - =M r -M r , = Y, {Yk+i - W k+1 ) 

k=r' 

for any r > r' > p. 



We can now state and prove a dual representation for the general multiple stopping problem in 
terms of martingales. 

Theorem 3 (Dual representation). In the setting described above, we have that: 
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(i) For any < i < d and any set of martingales [M r ' fc_1 ) , where 1 < k < L, 

^ ' r>j k -i 

and i =: jo < j\ < ■ ■ ■ < jjfc_i, ^ ZioWs 

y " * E * k^?^ (** - + E - 1 - . ( 3) 



fc=i 



fu,) ZioWs /or i > 



Y* L = max 

i<ji<—<3L<d 



v. fe=l / 



w/iere /or 1 < k < L, and i =: Jq < j\ < ■ • • < Jfc-i, \M* L fc+lj ' 1 '—' J ' fc 1 j j S a £) 00 & 

\ ' r>j k -i 

martingale of ( Y ; L - k+1 ' jl '~> jk -A 

\ ' r>j k -i 

Proof, (i) For the martingale family as stated we have for any chain of stopping times < r 1 < 
• • • < t l < d, 

L 

k=l 
L 



fc=l 

E EjE T fc_i (M^^- 1 -M^ 1 ^" 1 ) =0, 



fc=i 



hence, 



L 

*i _ c.,-,^ TIT I Y A \^ / fl / fi-fc+l>T 1 v,T" fc ~ 1 A# i_fc + 1 '' rl '--->' rfe 



= snp Ej X t1 _ tL + ^ M r ^ +1 ' r - M 

i< T i<...< T -t 



\ fc=l 

from which (i) follows directly. 

(ii) For any chain i < j± < ■ ■ ■ < jx, < 9 we may write (recalling jo := i) 

„ + E «_; fc+1J] h : - < fc+1J] h ) 

fe=i 

L 

V i f v *L-k+l,j 1 ,...,j k _ 1 v *L-k+l,ji,...,j k -i 



k=i 

l jfc-i 



fe=l l=3k-l 
L 

_ y*L i \ " / y*L-k,j 1 ,...,j k _ y*L-k+l,ji,...,jk-i 
fe=l 

+ E E (E,!^"*^ 1 '-"^- 1 - If £ " fc+1Jl '™ Jfc - 1 ) . (4) 
fe=i i=jk-i 
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By the reduction principle ^ it follows that (y r * L k+irfu— >i* i\ - g a SU permartingale 

which dominates the (virtual) cash-flow y* ' 1 '"' , ' Jfe ~ 1 ' T ' for k = 1, ...,L. Hence, expression Q 
is less than or equal to y* 1 '. It thus follows that 



max I X n JL 



L \ 

i f Tur*L-k+l,h,...,j k - 1 _ ^*i-fc+l,ji,...,i fe _ 1 \ \ < y *L 

\ ifc-i ife / J — * ' 

fc=i / 



i<jl<-<jL<9 

and then, an application of (i) finishes the proof. □ 

A straightforward consequence of Theorem [3] is the following dual representation in terms of 
approximate Snell envelopes. 

Corollary 4. For any set (Y r ^ ^~ r± ' J1, — ,Jk - 1 j j approximations to the Snell envelopes 

) . . with y ' J 

'r>j k -i 3L 



y; 1 - h - Lr - h ) ^ . with v n " " 



Y* L < Yt + Ei max > ; | V.'/ •" - Y 4 



i<jl<-<3L<d 



-'), 










■■■JL> 


L 




E 




fe=i 





A- 



J'fc-1 \ 
+ (E I y { j7 fc+1Jl, -" ifc - 1 - y^-Wi.-^-iA 



(5) 



Equality holds when the Snell envelopes are plugged in. 



Proof. Given (y L + 1 >Ji>— ) 3k i j we denote a corresponding family of Doob martingales 

^ ' r>j k -i 

by fMr^ 1 ^ 1 ""'^- 1 ) . Following the same manipulations as in @ and recalling that by 



definition Yj' 3 ' 1 '"' ,3L = Xj 1: ...,j L , we get 



fc=i \ 



L-k,j 1 ,...,j k _ yL-k+l,ji,...,j k -i 
3k 3 k 



_|_ ^ ^y-^-fc+ijli-Jfb-l _ yi-Hlji,.-jfc-i 



L 



/c=i 

Hence, the assertion is a mere reformulation of Theorem O 



i-fc+ljivjfc-i nA -L-k+l,j 1 ,...,j k _ 1 
3k 

k=l 



□ 



At this point, we stress that the dual representation from Theorem [3] relies on families of martin- 
gales (Mr +1J1 '- Jfe - 1 ) whose size is parametrized via the (k— l)-tuples . . . ,jk-i), k = 

^ ' r>j k -i 

1, . . . ,L. Hence, depending on the time horizon T and the number of exercise rights L a huge 
number of martingales M L ~ k+l,J1, - ,3k - 1 , k = 1, . . . , L, < j\ < . . . < < 9 = T + 1 is 
required in order to compute an upper price bound by the above dual formulation. It is thus of 
great importance to single out situations, in which a family of optimal martingales can be con- 
structed from a much smaller family of auxiliary processes. This will be the topic of Section [3j 
A motivating example in this respect is the standard multiple stopping problem. 
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Example 5 (Standard multiple stopping). Let Z be a nonnegative adapted process with 
Zj = for j = d, i.e. no penalty is imposed for unexercised rights. The standard multiple 
stopping problem is to maximize EEfc=i Z T k] over the set of ordered stopping times t 1 < • • • < 
t l such that T k < r k+1 or r k = r fc+1 = d. This means that at most one right can be exercised 
per day, but an arbitrary number of rights can be left unexercised, i.e. is exercised at time d. 
This problem can be put into our general setting by considering the cashflow 



ELi z ik > if ij+i = H ^ij = d. 
—N, else , 



for N € N. Note that the Snell envelope Y* L does not depend on the choice of N, because it is 
never optimal to exercise X in a way which gives a negative payment. Hence, letting N tend to 
oo , Theorem [3] yields, 

jk=jk+i^ jk=& k-1 / 

where for 1 < k < L, i < ji < ■ ■ ■ < jk—Xi \M* L fc+1,J ' 1 '*" ,J ' fc - 1 j i s a Doob martingale of 



r>jk- 

fc-1 L 



Vr J: n ■ = £ Z jp + sup E r £ Z 

p=l r<T k <-<T L p=k 

t p=tP+ 1 tP=8 and T k =jk-i ^ jk-i = d 

for r > jk~i- Define, for r > 0, 



■r 

y<7-fc...< 7 -£' 
t p =t p+i => t p=(9 




and denote the Doob martingale of (Y* L k+1 ) r >o by (M* L fc+1 ) r >o- -4s 



Y*L-k+l _ y*L-fc+l,ii,...j 



fc-i 



is J- j k l -measurable for r > jk-i, we can conclude by Remark^ that (M* )r>j k -i i> s a Doob 
martingale of (y* L fc+lj ' 1 '—' 5 ' fc - 1 j _ Hence we end up with the dual representation 

\ ' r>j k -i 

/ L L 

Y* L = max ( V Z ilr + V ( M* L ~ k+1 - M* L ~ k+1 
'■h<-<3L \ ^ ^ V Jk - 1 Jk 

;,. , , i,.=f) \k=l k=l 



*<Jl 
jk=jk+i^jk=9 



of Schoenmaken (20 id) . Here the potentially large family of optimal martingales (M* L ~ k+1 ^ 1, '"' : ' k - 1 ) ) 
k = 1, . . . , L, < j\, . . . < jk-i < d, collapses, in fact, to a family of L martingales, namely the 
Doob martingales ofY* k . 
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3 Generic cashflow with additive and multiplicative structure 



We now introduce a generic cashflow structure for which the dual representation simplifies in a 
similar way than for the standard multiple stopping problem in Example To this end let us 
consider for each k = 1, L and I = 1, L — 1 two adapted processes U k and V . We define a 
"pre-cashflow" 



k-l 



k=l 1=1 

which is assumed to satisfy Xj lt „ j L > —N for some (possibly large) N £ N. Concerning the 
processes U k and V 1 , we suppose that Uf is integrable for every k = 1, . . . , L and i = 0, . . . , d, 
and that Vf is strictly positive and bounded from above for every 1 = 1,..., L—l and i = 0, . . . , d. 
The multiple stopping problem which we have in mind is to optimally exercise this pre-cashflow 
under some constraints on the set of admissible stopping times, which we now formulate. We 
first define an adapted volume constraint process v with values in {1, L} such that vt is the 
maximum number of rights one may exercise at t, and such that vq = L. In order to formalize 
this constraint, we introduce for p > 1 the mapping £ p which acts on a non-decreasing p-tuple 
(ji,---,j P ) by 

£ P (ji, -Jp) ■= #{r ■■ 1 < r < p, j r = j p }. 

Hence, £ p denotes the number of rights exercised at j p in the non-decreasing chain < j\ < 
• • • < jp < 9. Obviously, an ordered chain of stopping times r 1 < • • • < r L satisfies the volume 
constraint if and only if £ p (r , r p ) < v t p for every p = 1, . . . , L. The second constraint, which 
we want to impose, is a refraction period which specifies the minimal waiting time between two 
exercises at different times. We admit random refraction periods, i.e. at each time i, < i < d, 
we fix a stopping time p l taking values in {i + 1, d}. If at least one right is exercised at time i, 
then the refraction period constraint imposes that the next right must either be exercised at the 
same time (if consistent with the volume constraint) or otherwise no earlier than p % . A standard 
case is p l = {i + 5) A d, where 1 < 5 < T is deterministic. Both constraints can be summarized 
by the binary Fj p -measurable random variable 



Cp(ji , jp) •- 



1, Vi<k p : Siijt, ...,ji) < v h and Vi<k p : ji > ji-i ==> ji > p jl - 1 
0, else, 



which is equal to 1, if and only if the constraints are satisfied when exercising at the p times 
ji < ■ ■ ■ < j P - 

The dynamic multiple stopping problem which we now study is 



Y* L = sup Ej 

i<T 1 <-<T L <d 

C l {t\...,t l )=1 



k-l 



£^n^> 



.k=l 



(6) 



i.e. the supremum is taken over all stopping times with values in {«,... ,T, d} which satisfy 
the volume constraint and the refraction period constraint. This problem fits in our general 
(unconstrained) setting by considering the cashflow 




jl> if ClUi,,-, jh) = 1, 



X h ,..., jL = { „ else (7) 
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To illustrate our motivation for studying the previous cashflow, let us have a look at the following 
examples. 



Example 6 (Swing options). We extend the situation in Example [5] by imposing volume 
constraints and refraction periods as decribed above. Hence, we have 

V}:=1, Z = 1,...,L-1, j = 0,...,d, 
Uf :=Zj p = l,...,L, j = 0,...,8, 

where we recall that Z is a nonnegative adapted process with Zq = 0. The multiple stopping 
problem then becomes 



sup E 

r 1 <-<r 1 ' 
C l (t\...,t l )=1 



leading to 



£2 r 

,k=l 



X 



Hk=l Z 3ki ifCdjln-jL) = 1, 



31, — ,JL 



-N, 



else. 



Here any N € N can be chosen because Z is no nnega tive. A dual approach fo r this m ultiple 
stopping problem was studied by \Bender \2011b\ ) and \Aleksandrov and Hamblil If 20 id ) under 
volume constraints, but with unit refr action period, i. e. p 1 = i + 1. The case with non-trivial 
constant refraction period is treated in \Bende\ hOllA) . but only under unit volume constraint, 
i.e. Vi = 1. A typical problem in the context of electricity markets which leads to this type of 
multiple stopping problem is the pricing of Swing option contracts, in which volume constraints 
and refraction periods are often imposed. This option pricing problem will be explained in more 
detail in our numerical study in Section^ 



Example 7 (Exponential utility). Under the assumptions of the previous example we can 
also maximize the exponential utility of exercising the cashflow Z% L-times while obeying the 
constraints. Given the risk aversion parameter a 6 (0, oo) the corresponding multiple stopping 
problem becomes 



sup E 

T 1 <-<T i 

C l (t\...,t l ) = 1 

This problem fits in our setting by considering 

k-l 



X 



3i,—,3L 



Eti^IK" ifCL(h,,..,3L) = h 
i=i 

—N, else 



with 



V 



e~ aZ i > and ll] 



0, ifk = l,...,L-l, 

-e~ aZ \ ifk = L, 



for j = 0,1,..., d and N > 2. 
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Example 8 (Portfolio liquidation). Suppose a (large) investor on a illiquid market wants 
to sell out (liquidate) L shares of a stock during the period {0, ...,T}. We assume that Sj > 0, 
j = 0, ... ,T, is the virtual stock price process reflecting the stock price evolution in the absence 
of the large investor's trading. In the spirit of \Schied and Slvnkd kOli ). Section 3.1, we model 
the price impact of the large investor by a resilience function G which we here apply to the 



3k 

k — 1 shares were already sold at dates < j± < • 



-price. Hence, the log-stock price lnS^'"' ,Jk 1 at time jk of the sale of the kth share, where 

< jk-i, is given by 

fc-i 



lnS£'~ J *- 1 



lnS j k -J2 G ^ k -ft- 



1=1 



We here choose the capped linear resilience function G{t) = 6(1 — at) + for constants a,b > 0. 
Assuming a short time horizon T < 1/a , the investor is thus faced with a multiple stopping 
problem 

r L 



sup E 

T 1 <-<T L 



k=l 



which fits in our framework by applying, for < j\ < ■ ■ ■ < ji, < T, the cashflow 



k-l 



;3L 



= E ^r^ 1 = E 4 exp - E b (1 - a (,, 

k=l k=l V 1=1 

L k-l 

= e ^ exp t & (° j ' fc ~ ^ ( k ~ i )i n exp (~ a ^) 
/c=i 

=E^n^ 



jo) 



i=i 



k=l 



1=1 



with 



U 



Sj exp [b (aj — 1) (k — 1)] and V- := exp (—abj) . 



(Note that the cemetery time d is irrelevant in this setting and we can e.g. set Ug = 0, Vg = 1 
to make sure that it is never optimal to exercise at this time). 



Similarly to the situation in Example [5l we now introduce a family of auxiliary multiple stopping 
problems Y* L ~ k+1 , which are not parameterized by the times ji, ■ ■ ■ ,jk—i, & t which the first 
rights were exercised. We will then show that a family of optimal martingales for the original 
multiple stopping problem (|6|) can be constructed via the Doob decomposition of the auxiliary 
problems. This then leads to a simplified dual representation for (0), which can be implemented 
in practice even when the maturity T and the number of rights L are large. 
Define 

Y r * L - k+1 := sup |E ^11^4 ] ^ 

T k ,-,T L \p=k l = k I 

r<T k <-<r L and C L _ k+1 (r k ,...,t l )=1 

with the convention 5^*° := 0. The following proposition states the Bellman principle for this 
multiple stopping problem. 
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Proposition 9 (Dynamic program). For r > and 1 < k < L we have, 



fc+n— 1 p—1 fc+n— 1 

+1 



Y* L ~ k+1 = max M r Y*+r k+1 , max V E^TIVr + ' ^X^* L ~^ n 

I r+i l<n<tv A (L-fc+1) ' ^— ' - L - L P 



p=fc Z=fc Z=fc 



Proof. From ([8]) we derive straightforwardly, 



Y*L— fc+1 



max sup E r V C^TTT^f + V t#, ' " V q 

0<n<<v A (L-fc+1) r<T fe+n < ... <r i ^ 11 1 



fc+n— 1 p—1 L fc+n— 1 p—1 



C i - fe+ i(r,...,r,r fe +",...,r i )=l 



p=fc Z=fc p=k+n l=k l=k+n 



-k + l 

(/k+n—1 p—1 
E r y;+r fc+1 > max y c/ r p TT^ 
7 + 1 l<n<u r A (L— fc+1) I ^ f-f 

\ p=K / = fc 

fc+n— 1 L p—1 

+ n vj sup E r ( ^ v*, n ^ 

Z=fc p<~<T fe + n <-<T £ p =k+n l=k+n 

C i _ fe „ n+1 (r fe +' l ,...,r L ) = l 



A standard argument shows that 



L p-l 



sup E r U t? ' [ V r' 

pr< T X + n<...< T L \ =k+n l = k + n 

C L _ k _ n + 1 (T k + n ,...,T L ) = l 

L p-l 

= E r sup E„r( £ U*, ] r =E r y; f i - fc -"+ 1 , (9) 

pr < r fc + n<...< r L V p=fc+n Z = fc+n 

Ci,_ fc _„ + i(r fe +™,...,r L ) = l 

which concludes the proof. □ 

We now establish a crucial relationship between the Snell envelopes y* i_fe + 1 >ii!— Jfc-i an d y* i ^ fc + 1 
defined in (JBJ). The following Proposition shows that y* Z/ ~ fc + 1 >-? 1 '---Jfc-i ) parameterized by the 
jfc's, can be represented in terms of Y* L ~ k+1 which avoids the j k 's. Notice that, for k = 1, both 
Snell envelopes coincide by definition. 

Proposition 10. Suppose 1 < k < L + 1. Under the condition Cj~-i(ji, ...,jk-i) = 1, we /iai>e 
(%) /or r > jfc-i it holds 



=1 Z=l !=1 



(jjj Further it holds 



fe-i p-i 
v *L-fc+l,ji,...,i fe _ 1 _ ^ rrP TTt/' 
J'fc-i i P 11 ii 

p=l Z=l 

fc-1 (fc-l+n p-l fc-l+n 

(=1 I p=k l=k l=k 
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where the maximum runs over the J-j kl -measurable set 

N(j u . . . , j fc _l) := W, 0<n< {v jk _ x - £ k -i{ji, jfc-l)) A (L - k + 1)}. 

Proof. For k = L + 1 both assertions are implied by the conventions Y* ^ 1 ''"'^ = Xj li4 __j L and 
Y*° = 0. Hence we assume for the remainder of the proof that 1 < k < L. 

(i) Under ...,j k -i) = 1, we have for r > j fe _i 

y*L-k+l,j 1 ,...,j k _ 1 
j p— \ k—1 I L p 1 \ 

= EKU v l + U v l , sup Er ■ ( 12 ) 

p=l i=l i=l r<rk<-< T L y p=k l=k J 

CL{jl,---,jk-l,T k ,...,T L ) = 1 

As r > , we obtain, thanks to ([9]), 

L p-l 



SUP K\Y, UP T»U. V T 

r<T k <■ ■ -<r^ 

Cl()lr"A-l,T t ,...,T i ) = i 



p=fc Z = fc 



W*-4 . sup Er XX'IIv; 

^*-i< T *<...< r i, y p=fc z=fc 

CL_ fe+ i(rfc,...,r-L) = l 



p-l 

I 



r < r fc<...< r L i z=fc 

C i _ fc+1 (r fe ,...,r i ) = l 

= i { ,<^- 1} E ^ p t:" +1 + i { ,>p^- 1} ^ L " fc+1 - (is) 

Hence, by combining (fTZj) and (fT5|) we get (i). 

(ii) Given that the first (k — 1) rights have been exercised at times j± < ■ ■ ■ < the number of 
the remaining {L — k + 1) rights which are also exercised at time jk-i must be chosen from the 
-7^ ^-measurable set N (h,- ■ ■ Jk-i) = {n; < n < (v jk ^ - £ k -i(ji, jk-i)) A (L - k + 1)}. 
These are the only choices which obey the volume constraint at time jk-i- Hence, 

n£N(3i,---,3k-i) p Jk-i <T n+k < ... <T L 

C L ^ k ^ n ~ 1 (r k +",..~r L )=l 

where the time index jk-i appears (n + 1) times in the nth term. It then follows, for fixed 
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n G N{ji,. . .,jk-i), 



sup EX, 



p Jfc_l< r n + fc<...< T L 



Jli— Jfe-li'" iJfe-l)T 

C i _ fc _ n+1 (r fc +",...,r i ) = l 

fc— 1 p— 1 ft— 1 ft— 1+n p— 1 

p=l 1=1 1=1 p=k l=k 

k-1 k-l+n L p— 1 

+ik n v l, . -p e ^ n ^ 

«=1 Z=fc p 3 fc-i<r' I+fc <-<r L p=fc+n Z=fc+n 

C L -fc- n +l(r fc +™,...,r i ) = l 

fc— 1 p— 1 fc— 1 fc— 1+n p— 1 fc— 1 fc— 1+n 

= ixik+ik e "jlikl+ik n vl^x m > 

p=l i=l i=l p=fc «=fc Z=l Z=fc 

making again use of (|9|). This implies (fTTj) . 

□ 



3.1 Dual representation based on Doob decompositions 

The goal of this subsection is to prove and discuss the following simplified version of the dual 
representation from Theorem [3] for multiple stopping problems of the form (J6j) . 

Theorem 11. Suppose Y* is given by {6jj. Then: 

(i) For any set of martingales (M J f' _fc+1 ) r>0 , k = 1,...,L, and any set of integrable adapted 
processes (^4^~ fc+1 ) r>0 , k = 1, . . . , L, it holds for i > and with jo := i 



Y* L < Ei max [YufUv, 

ClUi,-Jl)=1 



L k-1 

I 

h 

k=l 1=1 



k=l 1=1 

(ii) For every i > it holds with jo := i 

/ L k-1 



L k-1 

, \ " TT i ( M L-k+l _ M L-k+l , to .L-fc+1 _ to . AL-k+1 



= max E^ITT^ 

Cl(ji,-Jl)=1 



k ft=l J=l 



L k-1 



+ S TT\V! ( M* L ~ k+1 - Mf~ k+1 + E,. A*f~ k+1 - E ? - , A*i~ k+1 

11 Jl \ 3k-i 3k 3k p>k-\ 3k-i pOk-i 

k=l 1=1 / 

where M* L ~ k+1 , A* L ~ k+1 are the martingale part and the predictable part of the Doob decom- 
position of the auxiliary Snell envelopes Y* L ~ k+1 in respectively. 

We here recall that the Doob decomposition of Y ~ +1 is the unique decomposition of the form 

Y*L— fc+l y*L— fc+1 _|_ j^/j-*L—k+l j^*L—k+l 

where the martingale M* L ~ k+1 and the predictable process A* L ~ k+1 start in zero at time zero. 
In order to prove Theorem [TT] we need the following auxiliary result. 
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Proposition 12. Under the assumption of Theorem \lll a Doob martingale ofE r Y*^_ k+1 , say 



lk-1 



M* L fc+1 , is determined for r > jk-i by 



Vr 



M* L - h+1 - M* L ~ k+1 = M; L ~ k+1 - Mf~ k+1 + E h , A *^ k+1 - E r A*i~ k+1 . 

r 3k-l r J k _ 1 3k-l plk-l r p>k-l 



Proof. Using the Doob decomposition we may write 
EY *L-k+l 

p>k-\\j r 

_ -n w \r*L—k+l , -ti v*L— fc+1 

= L. . ^ n f + M* L ~ k+1 - Mf~ k+1 - E r A*i~ k+1 + Af- fc+1 

Ufc-l<T"<//*- 1 } \ Jfc-i r 3k-i 1 p J k-i 3k-i 

i 1 (\r*L—k+l I i/*L-fc+l t\ / t*L— fc+1 /i*L— fc+1 I a*L— fc+1 

~y^L — fc+ 1 _j_ -\/r*L— fc+1 -\/r*L— fc+1 

-ti TP — fc+1 -n 4*L — fc + 1 | A*L — fc+1 

1 {j*- 1 <r<p , *-4^" ^fc-i 1 {r>p»*- 1 }^ 

_ T^*L— fc+1 , /1*L— fc+1 -n / A*L— fc+1 /!*£— fc+l\ /I 4 N 

M *L-k+l _ M *L-k+l _ E A *f-k+l_ / 15 n 

Since line (|14p is the sum of a J-fc_i-measurable random variable and a predictable process and 
line (|15p is a martingale, the proposition follows. 

□ 

We now can prove the dual representation. 

Proof of Theorem\ll[ (i) Suppose that, for k = 1, . . . , L, M L ~ k+1 is a martingale and A L ~ k+1 
is an adapted and integrable process. Then, the process M r ' 3u '" ' k ~ 1 defined for r > jk-i 



via 

fc-l 

.V//- j > : := TT V\ ( M r L ~ k+1 - Mf~ k+1 + E,- , A L ' k+1 - E r A L - k+1 

11 3l \ r Jk-l JK- 1 p J k-l ' p>k-\ 

1=1 

is a martingale due to the boundedness of the V l, s. By Theorem we have 

fc=i 



Ei ^- nli ' X A '/: J i 

%<3i<—<3l<o V 

L fc-l 



Z— ✓ 11 Jl \ Jfc-l 3k J k p J k-i 3k-l 

fc=l 1=1 



Ik-l-^^k-l 



with X as defined in (|7|) for sufficiently large N € N. Letting AT tend to infinity, we observe 
that maximization only takes place over those j\ < ■ ■ ■ < ji which satisfy Cl(j'i, . . . = 1. 
Plugging in the definition of X for those j\ < • • • < yields the assertion. 
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(ii) We now apply Theorem \3^(ii) for X as denned in ((7]) with sufficiently large JVgN. Letting 
N tend to infinity again and substituting the definition of X, we obtain, 

/ L k-1 L \ 



Ci(jl,...jL) = l Vfc=1 1=1 k=1 



Jk 



whenever (m* L h+1 >^'->^-A are Doob mar tingales of (y* L k+1 '^'-'^-A . By 

^ ' r>j k -i V / r>j k _ 1 

Proposition 1101- f z ) and Proposition [12] we can take 

j^.*L-k+l,j 1 ,...,j k _ 1 _ j^.*L-k+l,ji,... s jk-i 

jk-1 jk 
k-1 

= Y\Vj ( M* L ~ k+1 - Mf~ k+1 + Ef, A*f-- fc+1 - E,-, , J 4*^- fc+lN ) . 



1=1 



□ 



Theorem 1111 gives a straightforward generic way to calculate upper bounds for multiple stopping 
problems of the form §6§ at time i = via Monte Carlo by performing the following steps in a 
Markovian setting: 

1. Solve the dynamic program in Proposition for the auxiliary problems Y* L ~ k+l approxi- 
mately, and let Y L ~ k+1 , k = 1, . . . , L, denote the respective approximations. 

2. Perform the Doob decomposition of Y L ~ k+1 , k = 1, . . . , L, num e rically , e.g. by one layer 



of nested Monte Carlo as suggested by lAndersen and Broadid (|2004l ) in the context of 
options with a single early exercise right. 

3. Plug the processes which stem from the numerical Doob decomposition into the formula 
of Theorem [TT1- ( i) and replace the outer expectation by the sample mean. 

This program will be carried out in more detail in Section 2] in the context of Swing options. 
Notice that for a large maturity and a large number of exercise rights, the pathwise maximum in 
the dual representation of Theorem 1 1 1 1 runs over a huge set. We will now show that, due to the 
special structure of the payoff in ([6|), this maximum can be computed efficiently by a recursion 
over the time steps and exercise levels. 

Given any L-tuple of martingales M = (M 1 ,...,M ) and any L-tuple of adapted processes 
A = (A 1 , ... , A L ), define, for n = 0, . . . , L and i = 0, . . . , d, 

L~n /k-1 \ 

9™' L {M,A) := max^ J~] ( ff Vj+ n ) (u^ k - (Mf- n ~ k+1 - Aff^ _fc+1 



jo=i<h<-<3L- .. 

C L -n(h,-,jL-n) = l 




By Theorem llll 

y * L < e[9q' l (m,a)} 

for any pair of L-tuples (M, A), and 

Y* L = 9q' L (M*,A*) 



for an optimal pair of L-tuples (M*, A*). Generalizing a related formula in lBalder et al.l (|201ll ) in 
the context of flexible (or choi 
by the following proposition. 



the context of flexible (or chooser) caps, the expression 9®' (M, A) can be recursively calculated 
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Proposition 13. For every L-tuple of martingales M = (M 1 , . . . ,M L ) and L -tuple of adapted 
processes A = (A 1 , . . . , A L ) it holds for i = 0, . . . , T and n = 0, . . . , L, 



9™' L (M, A) = maxfe(M,yl)-(A^-M^), max T(Il V i 

+ ( f[ v t x+n ) (e n + u ' L (M, a) - {M L ~ n ~ v - M t L - n ~ u ) 



v /k-1 

l+n \ jj-n-\-k 



\X=1 / 
+A L-n-u _ E . A L-n~u 

with 

e%\M,A) = f^(f[vt n )us+ k . 

k=l \i=x J 

Proof. The formula for 6g' L (M,A) is obvious by definition. In order to prove the recursive 
formula, we denote 

L—n /k—1 \ 



k=i \i=i / 

in ( aL— k— n+1 in 1 aL— k— n+V 



Then, 

6^ L (M,A)= max { max F n > L (j , . . . , j L -n)\- (16) 

u=0,...,ViA(L-n) <jo=-=iv=Kji,+l<---<3L-n 

Ci- n (ji,...,ii-n)=i 



For v = we get 



■ • m ^ F n ' L (h, ■■■,3L-n) 

Ci-ii(jlr-jl-n) = l 



. max . F^Uo, . . . ,j L - n ) - (M£? - M^ n ) 

]0=l + l<]l<-<JL-n 
Cl— n.{jl, — ,3L— n)=l 

0^(M, A) - - Mi~ n ). (17) 
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For v > we obtain 

max F n ' L (j ,...,jL-„.; 

J0=-=>=«<> + l<-<?L-n 

Ci-nCji,... Jz,-n)=l 

i/ /fc-1 \ 

W+n I Tjn+k 



fc=i \/=i / 

+ (n i/ i i+n ) . , <( e i n ^: + " j k* - u i 

Vj.i /»-•■<> s>«s / 

!/ /fe-1 \ 

E II ^/ +n 

fe=l \Z=l / 

+ ( n v " +n ) {{A L p r u - n - ^iA L r~ n ) - (M L p r~ n - m % 



fe=i \i=i / 

/ v \ L~n / k—1 



L-n-k+X 
jk 



rL—v—n s 

pi i pi i \ pi 

a=i / 

f v \ L~n / k—1 



\\=i j i»= p ^"+i<- -<^-n fe=i/+1 Vi=^+i / v 

t "i-u-]/Ul'+ll"-Ji-ni=- 1 

i j / L-n-fc+l , ■< / <L- fc— n+1 it? AL-k-i 



L-n-fc+1 

jk 



v /k-1 \ 

^i+n I jjn+k 



L—v—n s 



fc=i \«=i / 

(v \ L-n-v /k-1 \ 

H-n-i/Ulr-JI-n-») = l 

i ji .r L— n— fc+1 i n / \L— k— n— v+1 it? AL—k—n—v+1* 

+M jk-i l{fc>i a ife>ife— 1}' pjfe— i - ^jfc-i^fc-! 

1/ /fe-1 \ 

e n ^/ +n ^ +fc 

fc=l \l=l ) 



fL—n—v—k+1 



i / 



\A=1 / 

Plugging this identity and (fT7|) into (fTUj) yields the assertion. □ 



3.2 Dual representation based on Snell envelopes 

In this subsection we present a simplified version of the dual representation in Corollary H] in 
terms of approximate Snell envelopes for the multiple stopping problem of the form It reads 
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as follows. 

Theorem 14. Suppose Y* is given by (G|) for some fixed 0<i<d. Let {Y k ) i<k<L be any set 
of integrable approximations to (X* fc )i<fc<i defined in (|8j). We then have, with the conventions 
j ■= -l, pio ■- i f an d Y° = 0, 

•y*L ^yL 

1 i 1 i 

l ( ife-i fc-i 

s E w<™su, e e n^K^ 1 -^—) 

fe— 1 k+n— 1 p— 1 



+w.> rK( 1£ „^ +1) { e ^rw 

fc+n-l "\ 

-I- TT E- v L ~ k ~ n+l \ — Y 
+ 11 ^fc^JfeV* J J 



fc+n-l 

Moreover, the righthand side becomes zero ifY k = Y* k , fork = l,...,L. 

Proof. Suppose < i < d is fixed and assume that integrable and adapted processes Y L ~ k+1 , 
k = 1, ...,L, are given which we consider as approximations of the Snell envelopes of the 
auxiliary multiple stopping problems Y* L ~ k+1 . Following the relationships for the Snell en- 
velopes Y* L ~ k+l and Y* L ~ >J*— i i n Proposition 1101 we define for k > 1 approximations to 

y*£-fc+lji,...j fe _i y j a 

r> ifc _ lj (18) 



p=l !=1 Z=l 



and (for r = j/^ 



y.£-fc+lji,...,j fc _i 

ifc-i 
fc-l p-l 

- e^ik 

p=i i=i 

fc-1 (fc-l+n p-l fc-l+n 

Clearly, we define, for fe = 1, Y L $ = Y L . 

Applying Corollary H] for X as defined in ([6]) and the above approximations we obtain, 

L ( 

r*L / vt , ip. ryr , av / \rL-k,ji,... ,jk xr L -k+hhv>3k-l 



Y *L < Y L + Ei max r y^W* _ yf- 



C L (h,-,jL) = l k=1 V 

J=Jfc-i / 

where we again observe that the pathwise maximum is attained on the set Cl(Ji, j'l) = 1 by 
letting N (in the definition of X) tend to infinity. 
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In order to prove the upper bound, it is, in view of (|20|) . sufficient to show that, for i < j\ < 
■ ■ ■ < 3l < 9 with ClUi, ■ ■, Jl) = 1 the following assertions are true: 
(i) If fe = 2, . . . , L and j k > jk-i or if A; = 1, then 

yL-fc,Jr'l,...,j fc _ yi-fc + l,jlvjfe-l 



fc— 1 / I k+n p— 1 fc+n 

11 31 \ 0<n<(vj, -l)A(L-k) ) Jfeii Jfc 11 Jfc Jfc P Jfe | J 



l=1 v 0<n< K .-l)A(L- fc ) u=fe ;=fc ;=fc 

(ii) If fe = 2, . . . ,L and j k = jk-i, then 

yL-k,jl,...,j k _ yi-*:+l,il,...jfe-l ^ Q 

jfc jk — 



L-k+1 

jk 



(iii) If fe = 1 and i < r < — 1, or if fe = 2, . . . , L and /O^* -1 < r < jfc — 1, then 



i=i 



(iv) For fe = 2, . . . , L and jfc_i < r < p 3k ^ 1 



E r Y^ 1 k+1 ' 3l, '"' 3k ^ 1 - Y r L - k+1 ' jl '-' 3k - 1 < 0. 



We first show (i). To this end suppose that fe > 2 and jk > jk—i- Then, £ k (ji, ■ ■ ■ , j k ) = 1, 
which implies N(ji, . . . ,jk) = {n; < n < (v k — 1) A (L — fe)}. Hence, by (TT91) , 



L-k,j!,...,j k 



fc-1 p— 1 fc— 1 | fc+n p— 1 fc+n 

Subtracting the defining equation (|18p for from the above expression, we obtain 

(i), because j k > p )k - 1 . For fe = 1, we again get £\{ji) = 1 and (i) follows in the same way, 



rL,0 

jl 

In order to derive (ii), we note that £ k {ji,. ■ ■ , jk) = £k-l{h, ■ ■ ■ Jk-l) + 1 for j k = jk-l- Thus, 



taking the definition Y-^ = Y£ into account. 



-yL—k,ji r ..,jk 
jk 

fc-1 p-1 

rl 
'jl 

= 1 1=1 



E^IK 



k—l ( k+n p— 1 k+n 

+1\VI max J V 17? TT^* + ' 1^%**^"' 

fe-1 p-1 



p=i i=i 



fc— 1 (k+n p—1 k+n 

+T\vj. max V^, TTW. t + " ~[Vj. M jk x Y L - k ~ n 

\} x 0<n<(i) jfe _ 1 -£' fc _ 1 (ji,...,j fc _ 1 ))A(L-fc+l) | ^ Jfc - ll - L Jfe ~ 1 11 Jfc 1 ^-1 



I p=fc l=k l=k 



yL-k+l,ji,...,j k -l _ yi-fc+l)il)-Jfc-l 

ife-i it 
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We next prove (iii). The case k = 1 is trivial in view of the definition of Y L >*. Hence, we assume 
that k>2 and p>k-t < r < j k - 1. Then, r + 1 > r > p^- 1 > j k _ x and, thus, by (JTSJ 

p=i z=i i=i 

p=l Z=l Z=l 

Taking the difference of both equations yields (iii). 

It remains to show (iv). For k > 2 and jfe-i < r < /9 3 * -1 , (|18j) implies 

E r y r i + - 1 H1J1 = ^t/?n^'+^^^ 

p=l z=l 1=1 
Finally, for k > 2 and r = by (fT51) and (fT9]l . 

fc-l p-i fc-l 

p=l Z=l Z=l 

fe-l p-1 

p=i i=i 

fe-l ( fe-l+n p— 1 fe-l+n 

„i-fe+l,ji,...,j fc _ 1 
~~ j'fe-i 



Hence, the asserted upper bound for YJ* L — Y/^ is shown. This upper bound is zero, if Y k = Y* k , 

y-*L 

Ok 

fc+n— 1 P— 1 fe+n— 1 



for fc = 1, L, because, by Proposition [9j y*^ k+1 is a supermartingale and Y* L k+l dominates 



— { e ^ik+ n ^%^ n+i }- 



l<ra<u.,- A(L-fc+l) 



p=k 

□ 



As a spin-off result from Theorem [TH we may write the following upper bound for Y*' L which 
avoids the computation of the recursive maximum from Proposition [13] (cf. Schoenmakersl 
(|20ld ) [Remark 3.3] for a related result in the context of the standard multiple stopping problem) . 

Corollary 15. Suppose all assumptions and all conventions of Theorem^Tj^ are in force. Then, 

V*L y-L 



I r=i ~ ^ ' fe=l 



(fc+n— 1 p— 1 fe+n— 1 
V C/fTTW TT VjE j Y%- k - n+1 I -if"' 
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where 

k 

Vmax : =I]™ Xy /- 

i>i J 

1=1 

Moreover, the righthand side becomes zero ifY k = Y* k , for k = 1, ...,L. 

Proof. It is straightforward to check that the upper bound in this corollary is actually an upper 
bound to the righthand side of the estimate in Theorem 1141 That the bound is still tight, i.e. 
that the righthand side becomes zero, if Y k = Y* k , for k = 1,...,L, follows from the same 
argument as at the end of the proof of Theorem [T5J □ 



4 A numerical example 

We provide a numerical example for the dual representation of multiple stopping problems in 
the context of swing option pricing. Throughout this section, we assume i = 0, i.e. we provide 
confidence bounds for the swing opti on price at time 0. Precise l y, we consi der a stylized sw ing 
option, similar to those considered in Meinshausen and Hamblv ( 20041 ) and Bender ( 2011al ). In 



our setting, the holder of a swing option has the right to buy a certain quantity of electricity in 
the period from j = 0, . . . , T, for a fixed strike price K > 0, subject to the restriction that the 
option allows up to L > 1 exercise opportunities under the volume constraints Vj, and where a 
refraction period has to be taken into account. Here we choose T = 50 and recall that d := T+l. 
The price of electricity, (St)t=o,...,T, is modeled by the following discretized exponential Gaussian 
Ornstein-Uhlenbeck process 

log(Sj) = (1 - ^(logCSj-x) -n)+fi + aej, S = s > 0, (21) 

where (ej)j=i t 1S a family of independent standard normal random variables and the param- 
eters are specified by 

a = 0.5, k = 0.9, fi = 0, s = 1. 

We set Sq = 0, which means that no penalty is imposed, if the holder of the option does not 
exercise all rights. The payoff of the swing option is then given by X in ([6]) with 

Vj:=l, l = l,...,L-l, j = 0,...,8, 

U p j :=Z j :=Z(S,-) := (Sj - K)+, j = 0,...,8, p = l,...,L. 

In our numerical study we assume that the strike price is K = 1. As volume constraints we 
consider the situation of a unit volume constraint Vi = 1 for i = 0, . . . , T and the situation of 
an off-peak swing option with V{ = 1 on weekdays and = 2 on Saturdays and Sundays. The 
refraction period which we impose is a constant refraction period, i.e. p % = (i + 5) A d for various 
choices of the constant 5 £ N. 

In this Markovian framework, we produce confidence intervals for the price of the swing option 
at time i = by applying the following steps. The procedure below can easily be generalized 
to the generic cashflow structure of Section 3, provided the problem has a Markovian structure. 
(For notational convenience we only spell out the algorithm for the swing option case.) 
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4.1 Implementation 

Step 1: Precompute an approximation of the continuation values. We employ least squares 
Monte Carlo regression to obtain an approximation to the continuation values 



Cf\S 3 ) 


:= E 


v*i 

*j+l 




= E 


v *i 


Sj 


C* s '\Sj) 


:= E 


v *l 




= E 


v*i 


Sj 



1*1, 



0. 



C* T S >\S T ) = 0, 

with I = 1, . . . , L, where here and in the following j+1 and j+6 are to be understood as j + 1 A d 
and j + 5 A d. Recall that (^*')j=o,...,T is given by the dynamic program from Proposition [9j 
We simulate JVi independent paths (S™)™!^''''^ 1 ■ Choosing as basis functions 

ipi{x):=x, i>2{x) := (x - K) + , 

we use {S^Y^Zq'"'^ 1 in a straightforward least squares regression procedure to solve the dynamic 
program approximately, replacing the conditional expectations by the least squares Monte Carlo 
estimator. This yields approximations to C* 1,l (-) and C* S,l (-), denoted by Cj'\-) and C 5 -'\-). 

Step 2: Compute lower bounds. Given the functions C^'\-) and Cj'\-), we define a (suboptimal) 
stopping rule ( r J'')i<p<t ^ or ^ — 3 — ^ a l° n g a given trajectory (Sj)./=o...,T (which we suppress 
in the notation below) using the following iteration. Here t p ' 1 is interpreted as the time at which 
the investor exercises the pth right, if / rights are left at time j. 

Tj :=J-S; 
p:=k:= 0; 
while (p < I) do 



p+l,l 



:= inf { (rf + 5)Ad<r<8: max ( nZ r + C*' 1 -*™) > C}.' l ~ p }, 

I 3 l<n<v r A{l-p) V / 



p+l,l 
S:=T J > 



k := argmax!^^^^^) ( nZ s + C 



i8,l—p—n 



p+l,l p+2,l 

p := p + k, 



p+k,l 



(22) 



end 



When Cj' l (-) and Cj ,l (-) are replaced by C* 1,; (-) and C* S '\-), then this family of stopping times 

is optimal. Hence, ( r J' i ) 1 <p<f is a good family of stopping times, if the approximations of the 
continuation values in Step 1 are reasonably close to the true continuation values. 



Remark 16. (i) In the situation of unit volume constraint (i.e. v = 1), the stopping rule 
simplifies to r?' = j — 8 and 

rf = inf {{t p ~ 1 ' 1 + 5)/\d<r<d:Z r + Cf.' l ~ p > C}: l ~ p+l }, 1<1<L, l<p<l, 
compare with Eq. (3.7) in lBendeA koild ). 



(ii) In the situation of a trivial refr action period (i.e. 5 = l) l the ab ove construction of approx- 
imate stopping rules is also used in XAleksandrov and Hambh hoi A ). 
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Setting 

i i i 



p=l p=l p=l 



we have that Y} is a lower bound for Y * 1 . By the tower property of the conditional expectation, 
we also have 

l i 

p=l p=l 

As for simulations, we generate a new set of N2 independent paths of the underlying price process, 
which we again denote, in abuse of notation, by (SJ l )J i z r l, "'^ f2 . Along theses N2 trajectories we 

compute Tq'' and apply the notation 

T f' m , l<p<l, \<m<N 2 . 

Now the lower biased estimate Y_ for Y is calculated by averaging over the N2 realizations of 

1 N * 1 

^=F^^^ a 1</<L. (23) 



<2 

m=l p=l 

Similarly, we also construct approximations 

N 2 l , N 2 I 



E, 



m=l p=l m=l p=l 

of E0Z1 and E rJ>, which we store for later use. 

For constructing confidence intervals, we also save the empirical standard deviation stddev(Y ). 

Step 3: Compute approximations to the Snell envelopes. Using the stopping rule ()22|) . we 
consider a family of random variables 

I 

Y< := Ej Z rP , h 1 < I < L, < j < T, (24) 



p=i J 



/ ,\i<i<i 

which is an approximation to the Snell envelope I Y* . We apply the following procedure 



Jo<j<T 

to simulate Y^-: 

We simulate a new set of A3 paths of the underlying {SJ 1 )^™^ 3 (abusing the notation, again). 
We refer to these paths as the outer paths. We now fix a pair [m, j) and compute approximations 

of Yjj, EjYj +1 , and EjY_^ +( j along the mth outer path which are denoted by Y_- , W^Y}^ +1 , and 

Kj^Yj+g, respectively. In these approximations the conditional expectations are replaced by 
the sample mean over a set of inner simulations. Hence, for the fixed path S m and the fixed 
time point j, we generate independent sample paths of (SV)r=j,...,T under the conditional 
law given that S,- = <S"\ These inner paths are denoted by (S^Yz]''"'^ A , suppressing here and 
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in the following the dependence on (m, j). Along the inner paths S u we compute the stopping 
times rf' for i = j, j + l,j + S in (|2"2"j) and apply the notation 



rf 1 < p < I, 1 < I < L, u = l,...,N 4 . 



We now define 



l N 4 l 

Similarly, we approximate EjY^ +1 for the fixed j along the fixed mth outer path by 

N 4 i 



iV 4 — — ■ ' J+ i 

z/=l p=l 

taking the tower property of the conditional expectation into account. The approximation 
E™Y^- +(5 is obtained analogously. 

Remark 17. Note that, for j = approximations Y_ , EqY\, and EqY^ °f — d' ^oHi> anc ^ 
EoH^ were already obtained based on the N2-samples in Step 2. As typically N2 > N4 these 
approximations are more accurate. Hence, one can perform Step 3 for j > 1 only and set 

£o m :=£o, E^-EoFi, Eg*Y$ := E F<, m = l,...,M. 

TTiis iric/c 0/ applying the more accurate non-nested Monte Carlo simulation of Step 2 at time 
leads to a significant decrease of the variance in the simulation of the upper bound. This is 
in the same spirit as the computation of l ow variance upper bounds for the standard stopping 
problem from lAndersen and Broadiei (200$)). 



Step 4: Compute the upper bounds. The Doob decomposition of Yj yields the pair (M^A^). 
Note that due to 

Ml + x-Ml=Y! i+1 -E i Y! i+1 , 

and 

-(M! i+5 - Mi) + M+s - M+« = EiZU - rU 

we can rewrite the recursion formula in Proposition [13] as 

nn,L 



max 



Uti + EtYf-f - Yjtf, max (yZ { + 6^ L + E^T" " } 

l<i'<tiiA(L— n) \ / J 



We now introduce approximations Q™> L > m Q f along the mth outer path of Step 3 by replac- 
ing Y*., EjY l j+1 , and EjYj j+s with their simulated counterparts Y_- m , E™Y^ +1 , and EJY}- +5 
constructed in Step 3. 

As simulation based estimate for the upper bound, we use 

-tritp,L \ r /jO,L,m 

' _ "/Vo • 



N 3 

J m=l 
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Replacing the conditional expectations by the sample mean in 6 ' ,m introduces an additional 
bias up thanks to Jensen's inequality and the convexity of the maximum. Hence, the estimator 
Yq P ' L is biased up by Theorem [11] and Proposition fT3l 

Finally, a 95% confidence interval on the price of the swing option is given by 
Yq - 1.96 x stddev(£o),y o up ' L + 1.96 x stddev(Y up ' L ) . 











95% confidence 








95% confidence 


5 


L 






interval 


L 


ti 




interval 




1 


2 


3.3116 


3.3211 


[3.30738, 3.32229] 


3 


4.53627 


4.54806 


[4.53118, 


4.54938] 


2 


2 


3.27513 


3.28469 


[3.27094, 3.28587] 


3 


4.43753 


4.45154 


[4.43252, 


4.45295] 


3 


2 


3.2525 


3.26286 


[3.2483, 3.26414] 


3 


4.36706 


4.38245 


[4.36204, 


4.38392] 


4 


2 


3.2313 


3.24083 


[3.22716, 3.242] 


3 


4.29996 


4.31656 


[4.29502, 


4.31813] 


5 


2 


3.20906 


3.22061 


[3.20496, 3.22199] 


3 


4.29996 


4.31656 


[4.29502, 


4.31813] 


6 


2 


3.18613 


3.19809 


[3.18197, 3.19948] 


3 


4.15557 


4.17514 


[4.15063, 


4.17697] 


8 


2 


3.13625 


3.14984 


[3.13213, 3.15143] 


3 


3.99773 


4.01954 


[3.99289, 


4.02158] 


10 


2 


3.09022 


3.10332 


[3.08613, 3.1048] 


3 


3.83377 


3.8528 


[3.82898, 


3.85464] 


12 


2 


3.03874 


3.05196 


[3.03468, 3.05356] 


3 


3.65492 


3.67658 


[3.65023, 


3.67868] 


14 


2 


2.98727 


3.00048 


[2.98321, 3.00199] 


3 


3.47017 


3.49061 


[3.46558, 


3.49258] 


16 


2 


2.92751 


2.94214 


[2.9235, 2.9438] 


3 


3.27524 


3.29482 


[3.27077, 


3.29674] 


18 


2 


2.87368 


2.8888 


[2.86964, 2.89049] 


3 


3.09209 


3.11002 


[3.08775, 


3.11186] 


20 


2 


2.81521 


2.83005 


[2.81123, 2.83173] 


3 


2.91951 


2.93649 


[2.91536, 


2.9383] 



Table 1: Unit volume constraints (vj = 1). Numerical results based on the approximation Y- 
to the Snell envelope via the stopping rule (|22p for two and three exercise rights. 



4.2 Numerical results: swing options with unit volume constraints 



We now present some numerical results which the above algorithm produces for the swing 
option contract as specified at the beginning of this section. Let us first consider the situation 
of a unit volume constraint, i.e. Vj := 1 for j = 0, ...,T. We recall that 5 € N denotes a 
constant refrac tion per i od. In this setting the dual representation of Theorem [IT] reduces to the 
one derived in iBenderi (j2011al l. In the latter paper the same swing option example is treated 
numerically but for up to three exercise rights only. Thanks to the new recursion formula in 
Proposition [13] we can now efficiently treat the case of a large number of exercise rights (here up 
to L = 10). Moreover, t he upper boun d algorithm in IBenderi (|2011al ) differs slightly from the 
one we propose here. In Bender ( 2011al ) the upper bound is calculated based on the numerical 
Doob decomposition of 

Y^max^ + Cf- 1 ,^} 

while we here utilize the numerical Doob decomposition of Yj := Ej Ylp=i ^ t p,i- 
The choice of simulation parameters in our study is as follows: in Step 1, we choose N± = 1000 
paths for the least squares Monte Carlo regression to approximate the continuation function. 
In Step 2, the lower bound is simulated using iV~2 = 300000 paths and in Step 3, we employ 
N$ = 2000 outer and N± = 100 inner paths for the computation of the upper bound. Moreover, 
we use the variance reduction method from Remark 1171 
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Table [T] depicts the numerical results for the case of two and three exercise rights for a refraction 
period ranging from 1 to 20. We observe that the relative length of the 9 5%-confidence in tervals 
is less than 1% in all cases. A comparison with the numerical results in iBenderl l[2011a] ) shows 
that the differences in the upper price estimator based on Y l and Y_ l are negligible, but the 
variance reduction method of Remark [T7] shrinks the confidence interval significantly. 
The numerical results for the case of a larger number of exercise rights (L = 4, 6, 8, 10) are 
presented in Table [2j Due to the time horizon of 50 days, it may happen that, for a large 
number of rights and a large refraction period, some exercise rights cannot be used by the 
investor. This explains why e.g. the price bounds for the swing option with refraction period 
5 = 14 are the same for L = 4 and L = 6 rights. Concerning the accuracy of our numerical 
procedure we emphasize that the relative difference between lower and upper bound is still less 
than 1% even in the case of 10 exercise rights. 











95% confidence 








95% confidence 


5 


L 


2c 1 




interval 


L 


2f 


v up,L 


interval 




1 


4 


5.60136 


5.614 


[5.59554, 5.61527] 


6 


7.38677 


7.40107 


[7.37977, 


7.4023] 


2 


4 


5.41347 


5.43091 


[5.4078, 5.43249] 


6 


6.94554 


6.97364 


[6.93882, 


6.97562] 


3 


4 


5.27248 


5.29342 


[5.2668, 5.29509] 


6 


6.58739 


6.62054 


[6.58071, 


6.62272] 


4 


4 


5.13119 


5.15543 


[5.12562, 5.15741] 


6 


6.2151 


6.25165 


[6.2086, 


6.2541] 


5 


4 


4.98353 


5.01117 


[4.97802, 5.0133] 


6 


5.81445 


5.85591 


[5.80811, 


5.85862] 


6 


4 


4.82479 


4.85239 


[4.81928, 4.85454] 


6 


5.4079 


5.44248 


[5.40174, 


5.44491] 


8 


4 


4.49057 


4.51822 


[4.48525, 4.52041] 


6 


4.68469 


4.71574 


[4.67908, 


4.71808] 


10 


4 


4.13658 


4.16231 


[4.13141, 4.16444] 


6 


4.1662 


4.19164 


[4.16098, 


4.19373] 


12 


4 


3.78981 


3.81429 


[3.78491, 3.81652] 


6 


3.78992 


3.81448 


[3.78502, 


3.81671] 


14 


4 


3.50023 


3.52138 


[3.49558, 3.52334] 


6 


3.50023 


3.52138 


[3.49558, 


3.52334] 


1 


8 


8.83286 


8.84907 


[8.82488, 8.85034] 


10 


10.0219 


10.0391 


[10.0131, 


10.0404] 


2 


8 


8.04508 


8.08421 


[8.03754, 8.08651] 


10 


8.80264 


8.85093 


[8.79443, 


8.85353] 


3 


8 


7.36943 


7.41474 


[7.36205, 7.41734] 


10 


7.73096 


7.78117 


[7.72312, 


7.78394] 


4 


8 


6.66726 


6.71129 


[6.66021, 6.71389] 


10 


6.74649 


6.78596 


[6.73929, 


6.78847] 


5 


8 


5.99864 


6.0388 


[5.99202, 6.0414] 


10 


6.0035 


6.04305 


[5.99687, 


6.04558] 


6 


8 


5.45188 


5.48518 


[5.44563, 5.48748] 


10 


5.45187 


5.48518 


[5.44563, 


5.48748] 



Table 2: Unit volume constraints (vj = 1). Numerical results based on the approximation Y_ 
to the Snell envelope via the stopping rule (|22p for a higher number of exercise rights. 



4.3 Numerical results: off-peak swing option 

We now consider a swing option which allows for buying at most one package of electricity on 
weekdays and two packages on Saturdays and Sundays (off-peak period). Hence, we have for 
j = 0, . . . , 50 the volume constraints 

v _ = fl) if J is a week day, ^ 
[2, if j is a weekend day, 

where we start in j = on a Monday. 

We run the above algorithm with JVi = 10000, N 2 = 300000, N 3 = 2000, and iV 4 = 100 sample 
paths. The numerical results for this off-peak swing option are presented in Table for various 
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£ L 




95% confidence 








95% confidence 


<5 


L 


yUp,L 


interval 




L 


25 




interval 


1 


2 


3.39804 


3.40779 


[3.39342, 


3.409] 


3 


4.72241 


4.73543 


[4.71667, 4.73682] 


2 


2 


3.36368 


3.37544 


[3.35908, 


3.37676] 


3 


4.63568 


4.65304 


[4.62998, 4.65468] 


3 


2 


3.34728 


3.35873 


[3.34265, 


3.36017] 


3 


4.58105 


4.59793 


[4.57532, 4.59971] 


4 


2 


3.32763 


3.34072 


[3.32302, 


3.34223] 


3 


4.52367 


4.54475 


[4.51799, 4.5467] 


5 


2 


3.30829 


3.32101 


[3.30366, 


3.32246] 


3 


4.46437 


4.48607 


[4.45867, 4.48806] 


6 


2 


3.28626 


3.29915 


[3.28162, 


3.30067] 


3 


4.40272 


4.42552 


[4.39701, 4.42758] 


8 


2 


3.24383 


3.26104 


[3.23921, 


3.26284] 


3 


4.29259 


4.31753 


[4.2869, 4.31985] 


10 


2 


3.20872 


3.22244 


[3.20409, 


3.22404] 


3 


4.18401 


4.2091 


[4.1783, 4.21161] 


12 


2 


3.16994 


3.18626 


[3.16529, 


3.18807] 


3 


4.06657 


4.09375 


[4.06086, 4.09637] 


14 


2 


3.12722 


3.14467 


[3.12253, 


3.14666] 


3 


3.95351 


3.97621 


[3.94778, 3.97861] 


16 


2 


3.08634 


3.10272 


[3.08169, 


3.10464] 


3 


3.86402 


3.88597 


[3.85829, 3.88827] 


18 


2 


3.05186 


3.06793 


[3.04713, 


3.06991] 


3 


3.75608 


3.77801 


[3.75036, 3.78047] 


20 


2 


3.01557 


3.03044 


[3.01088, 


3.03233] 


3 


3.65138 


3.67244 


[3.64572, 3.67486] 


1 


4 


5.89744 


5.91312 


[5.89078, 


5.91464] 


6 


7.91351 


7.93364 


[7.90538, 7.93527] 


2 


4 


5.73736 


5.76003 


[5.73078, 


5.76192] 


6 


7.55394 


7.58547 


[7.54595, 7.58773] 


3 


4 


5.62688 


5.65097 


[5.62027, 


5.65305] 


6 


7.28486 


7.32349 


[7.27684, 7.32611] 


4 


4 


5.51763 


5.5468 


[5.51105, 


5.54915] 


6 


7.01995 


7.06275 


[7.01198, 7.06577] 


5 


4 


5.40154 


5.43295 


[5.39496, 


5.43546] 


6 


6.7418 


6.78724 


[6.73383,6.79042] 


6 


4 


5.27976 


5.30967 


[5.27316, 


5.31227] 


6 


6.45197 


6.49774 


[6.44401, 6.50102] 


8 


4 


5.06733 


5.10055 


[5.06078, 


5.10335] 


6 


5.9401 


5.98546 


[5.93238, 5.989] 


10 


4 


4.85039 


4.88637 


[4.84386, 


4.88953] 


6 


5.46672 


5.50782 


[5.45916, 5.51116] 


12 


4 


4.63227 


4.66583 


[4.62575, 


4.66884] 


6 


5.08473 


5.11779 


[5.07729, 5.12079] 


14 


4 


4.43104 


4.45997 


[4.42453, 


4.46279] 


6 


4.76734 


4.79957 


[4.76006, 4.80269] 


16 


4 


4.25079 


4.27955 


[4.24441, 


4.28237] 


6 


4.39537 


4.42311 


[4.38864, 4.42584] 


18 


4 


4.07804 


4.10338 


[4.07164, 


4.10605] 


6 


4.18139 


4.20744 


[4.17473, 4.21004] 


20 


4 


3.94562 


3.96789 


[3.93923, 


3.97034] 


6 


4.02465 


4.04716 


[4.01803, 4.04968] 


1 


8 


9.60253 


9.62348 


[9.59318, 


9.62507] 


10 


11.0436 


11.0661 


[11.0332, 11.0677] 


2 


8 


8.97188 


9.01806 


[8.96279, 


9.02078] 


10 


10.0822 


10.1411 


[10.0721, 10.1443] 


3 


8 


8.48335 


8.53629 


[8.47425, 


8.53952] 


10 


9.31393 


9.3793 


[9.30393, 9.38283] 


4 


8 


8.00789 


8.06203 


[7.99887, 


8.06551] 


10 


8.58082 


8.63832 


[8.57102, 8.64178] 


5 


8 


7.5251 


7.57926 


[7.5161, 


7.58278] 


10 


7.9058 


7.961 


[7.89611, 7.96454] 


6 


8 


7.06562 


7.11754 


[7.05669, 


7.12102] 


10 


7.33533 


7.38481 


[7.32577, 7.38835] 


8 


8 


6.18418 


6.23051 


[6.17596, 


6.23403] 


10 


6.20274 


6.24781 


[6.19445, 6.2513] 


10 


8 


5.54885 


5.58774 


[5.54107, 


5.59089] 


10 


5.54885 


5.58774 


[5.54107, 5.59089] 



Table 3: Off-peak volume constraints. Numerical results for the off-peak swing option for various 
exercise rights and refraction periods. The simulations are based on the approximation Y_- to 
the Snell envelope via the stopping rule ([22]) . 



26 



choices of the number L of exercise rights and the length 5 of the refraction period. Notice 
that the dual representations yet available in the literature do not cover the case of a nontrivial 
refraction period (5 ^ 1) in combination with nontrivial volume constraints (v ^ 1). Due to 
the feature of allowing for exercising twice on weekends, the swing option prices are now higher 
than in the example with unit volume constraint. Moreover, additional rights can now become 
beneficial in situations in which they could not be exercised under the unit volume constraint 
(e.g. the additional 8th right when the refraction period is 5 = 14). As for accuracy, we again 
observe that the relative length of the 95% confidence interval is less than 1% in all cases, which 
demonstrates that the algorithm performs equally well in the presence of volume constraints. 



<5 


L 


v up,L 


upper bound usine Bender ( 201 lb) 


1 


1 


1.86485 (0.0019) 


1.8638 (0.0019) 


1 


2 


3.40832 (0.003) 


3.4078 (0.003) 


1 


3 


4.73509 (0.0037) 


4.7368 (0.0038) 


1 


4 


5.90956 (0.0043) 


5.9170 (0.0045) 


1 


5 


6.96665 (0.00047) 


6.98 (0.0052) 


1 


6 


7.92669 (0.005) 


7.9470 (0.0058) 


1 


7 


8.80743 (0.0055) 


8.8327 (0.0062) 


1 


8 


9.61643 (0.0058) 


9.6493 (0.0069) 


1 


9 


10.3642 (0.0061) 


10.4040 (0.0074) 


1 


10 


11.0553 (0.00064) 


11.1035 (0.0079) 



Table 4: Off-peak volume constraints. A compar i son bet ween our upper bounds and the upper 
bounds obtained via the algorithm from Bender ( 2011bl ) for the case of unit refraction period. 
Standard deviations are displayed in parentheses. 



In the case of unit refraction period 6 = 1, upper p rice bounds for the off-peak swing option 
can also be computed by the dual representation of Bender ( 2011bl ) for the marginal price of 
a mul tiple exercise option. This approach generalizes the ideas of iMeinshausen and Hamblv 
(|2004l ) : An upper biased estimate for the marginal price of having an additional Ith right is 
computed in terms of one martingale and (I — 1) stopping times. By summing up these upper 
bounds for the marginal prices, one finally ends with an upper biased estimate for the option 
price. This approach is based on the fact that, roughly speaking, under the assumption of a 
trivial refraction period (5 = 1) optimal exercise times for the problem with (I — 1) rights are 
also optimal for the problem with / rights, if one adds one additional exercise time in a clever 
way. This is clearly not possible in general in the presence of a nontrivial refraction period. So it 
seems that this alternative approach cannot be easily generalized to include refraction periods. 
Table H] compares the upper bounds obtained using our method and the method from iBender 
(|2011bl ) for the unit refraction case 6 = 1. We mention that in Table HI the variance reduction 
method from Remark [T7] is not applied for both algorithms. As both methods are run with 
the same number of sample paths and the nested maximum in our method can be efficiently 
calculated by the recursion formula in Proposition [T3l the computational effort is roughly the 
same for both algorithms. We observe that, as the number of exercise rights increases, our 
method of directly tackling the Snell envelope pr oduces upper b ounds that become lower than 
the algorithm tackling the marginal values from iBenderl (|2011bl ). Whereas the differences for 
L = 1, ... ,4 are numerically not significant yet, they however become noticeable starting from 
L = 5 and are striking for e.g. L = 10. We also note that the larger L, the better our method 
performs concerning the variance of the upper bounds. At large, we conclude that if one is 
mainly interested in the price (and not the marginal price) of the swing option, our new method 
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performs better than the algorithm from iBenderl (|2011bh . Moreover, it is applicable to a larger 
class of problems. 
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